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Abstract 

We study instanton effects on B ^ Xg^f decay using the heavy quark ef- 
fective field theory and the operator product expansion. In the dilute gas 
approximation the effect is negligibly small. This is in contrast to the result 
of the instanton effect in inclusive semileptonic 6 — > n decay but similar to 
the inclusive hadronic r decay. We discuss the similarities and differences of 
the B decay compared to inclusive hadronic r decay and semileptonic 

B decay. 
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I. INTRODUCTION 



Rare decays of hadrons containing a single b quark are extensively studied as a sensitive 
probe to new physics beyond the standard model. At the electroweak scale the new physics 
introduces extra contribution to the coefficient functions of local operators in the low-energy 
effective theory ||T],01. Radiative 6 ^ 57 decay is of such type which is induced by penguin 
diagrams with virtual heavy quarks in the loop. Recent CLEO collaboration has reported 
the branching ratio of B ^ K*'~f as (4.5 ± 1.9 ± 0.9)% However theoretical estimates of 
the branching ratio are model dependent. So far all the theoretical estimates of the exclusive 
channels have been based on either the quark model or the QCD sum rules . The main 
uncertainty lies in evaluating the matrix element of the penguin operator between hadronic 
states. Therefore it is expected to be more reliable to consider the inclusive radiative B 
decays in which all the final states containing a single s quark are summed over. 

In calculating the inclusive decay rate we follow the approach of Chay et.al. to use 
the heavy quark effective field theory (HQEFT) and the operator product expansion (OPE). 
Bigi et.al. Neubert and Falk et.al. have performed the analysis and have found 
that the leading contribution is the parton-model result. They have systematically obtained 
corrections coming from the hadronic matrix elements of higher dimensional operators. The 
size of the corrections is about 5% of the parton-model result. There are other sources of 
corrections. Perturbative QCD corrections have been studied by Ali and Grueb |[T0|1 . 

One also expects additional nonperturbative QCD corrections. In the case of 5 ^ XuCVe 



decay we have calculated the instanton contribution using the dilute gas approximation ||Tl 



and have found that the instanton contribution to the total decay rate is large. In order to 
make the perturbation theory reliable, it is necessary to introduce a smearing prescription. 
In this paper we study the nonperturbative effects induced by instantons in B ^ Xg'j 
decay. We find that the instanton contribution to the total decay rate in this case is small. 
The striking difference of the instanton contribution in both decays may be surprising. We 
account for this difference by comparing the kinematics in the two cases. 
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This paper is organized as follows. In Section 2 using the OPE and HQEFT we formulate 
the decay rate from the effective Hamiltonian relevant to the B Xg'y decay. We relate 
the rate to the forward Compton scattering amplitude and discuss the analytic structure 
of the amplitude. In Section 3 we study the instanton contribution to the decay rate as a 
calculable example of nonperturbative corrections to the decay rate. We calculate explicitly 
the instanton contribution at leading order. In Section 4 we discuss the kinematics of the 
6 — > S7 decay and compare it to that of the b — > ueVe and the inclusive hadronic r decays. 
We find that the kinematics of the radiative S-meson decay is different from that of the 
semileptonic B decay but similar to that of the inclusive hadronic r decay. We explain how 
the kinematics affects the size of the instanton contribution in these cases. In Section 5 we 
summarize our results and give a conclusion. 

II. OPERATOR PRODUCT EXPANSION 

In the standard model the effective Hamiltonian for 6 — > s decay consists of various 
operators. They are generated after integrating out the W bosons, Higgs particle and heavy 
fermions at the weak scale. The effective Hamiltonian describing b ^ s decay can be written 
as 

= -^v,,v:,Y.<^,{^,)OA^^)■ (2-i) 

For inclusive b ^ sj decay only the operator C7 contributes, which is 

O7 = ^{m^LCJ^^R + msSRa^^bL)F^,. (2.2) 
The effective Hamiltonian for b ^ sj decay at = m{, is written as 

Hes = -^VtbV*cj{mb)07im,) (2.3) 

where the Wilson coefficient c-j{m},) is scaled down from ^ — Mw to /i — nif, using the 
renormalization group equation. 

The radiative S-meson decay rate is given by 



dT = ^\V,bV;f\crimb)\' ^ i2n)' 6^'\p b - Px - k)\{X,^\Or\B)\'diFS), (2.4) 

where the sum includes all the possible hadronic final states Xs containing an s quark and 
the photon polarizations A. = rribv'^ and are the four momenta of the B meson and 
the photon respectively. The phase space is written as 

diP^)= 77^.^ = T^Jdyy (2.5) 



(27r)3 2E^ 167r2 

where we have rescaled the photon energy as y = 2E^/mi,, < y < 1. Using the explicit 
form of the operator O7 and summing over photon polarizations we get 

dT = ^^\V,bV;f\cr{mb)\'ydy 
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X (2vr)^5(^) {pB -Px-k)\ (X„ 7| (m,Si + m,SR)a>^''F^,b\B)\^ (2.6) 

Xs,X 



The matrix element-squared in Eq. (|2.6|) can be factorized as 

Ex. {27ry6^'\pB - Px - k){B\ba^'''{mbSL + msSR)\Xs) 
X (X,|(mfeSL + m,s^j)a"^6|5) 

X ^(0|F^,|7;A,fc)(7;A,fc|F„^|0). (2.7) 

A 

The matrix elements of F^^ are written as (0|F^i,|7; k, A) = i^k^j^e^^^ (k) — k^e^^\k)) . Summing 
over photon polarizations, Y.x^\^'Kk)el^\k) = —g^u-, we have 

^(0|F^i.|7; fc, A)(7; fc, = {k^kag^^p - k^k^g^f} + k^^kpg^^ - Kkpg^^a)- (2.8) 

A 

The inclusive decay rate is related to the quantity W defined by 
W = {2'KfY.^^'\pB-px-k) 

Xs 

X (5|6[/;,7,](mbSL + m,Sij)|X,)(X,|(mbSL + m,SR)[^,7n&|5). (2.9) 

W depends on the photon momentum k^. In the complex y = 2v ■ k/rrib plane with fc^ fixed, 
W is related to the discontinuity of the forward Compton scattering amplitude T{y-k) across 
a physical cut as 



W = 2lmT. (2.10) 

For B decays it is appropriate to use the HQEFT in which the full QCD b field is 
expressed in terms of for 6-quark velocity v. The b^ field is defined by the heavy b field 
b^{x) = e*'^'^'^6(x) which satisfies '/)by = b^. We can apply the techniques of the OPE to 
expand T in terms of matrix elements of local operators involving b^ fields 

T{v ■k) = -ij d^xe'^'^'^-^^^BlTiJ^^ix) J^(0)}|5) 

= T.C:ik){B\Oi-^B), (2.11) 



n,v 



where 



= b4/c,-f^']{mbSL + msSR). (2.12) 



In Eq. ( p. 11 ) O^") are local operators involving b^by bilinears. The coefficient functions 
depend explicitly on k. Short-distance physics is contained in the coefficient functions while 
the matrix elements describe long-distance physics. 

In order to evaluate the decay rate it is necessary to examine the analytic structure of 
T. In the complex y plane the cut relevant to the decay is located on the real axis y < 1 —r 
where r = m^/m^. The cut extends down below y = 0, which is the boundary of the physical 
region for the decay process (in which the photon has a lightlike momentum). We will refer 
to the region along the real axis for 

0<y<l-r (2.13) 

as the "physical cut". There is also another cut 

y>(2 + v^)'-l (2.14) 

corresponding to other physical processes. These cuts in the complex plane are shown in 
Fig. 1. The discontinuity of T along this physical cut yields W for b ^ s'y decay. The 
physically allowed region for ?/ is < ?/ < 1 — r. We can introduce the cutoff of the photon's 
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energy such that He y ^ ^ — This cutoff is necessary experimentally since it is difficult 
to observe a soft photon. 

As ?/ ^ 1 — r, that is, at the right end of the physical cut, the invariant mass of the 
hadronic system approaches nts and we enter the region where resonances dominate. We 
expect large perturbative QCD corrections in this region. However we are able to reliably 
compute suitable averages of physical quantities by integrating around the physical cut 
along the contour C shown in Fig. 1. The contour integral along C and C contains no 
singularities, so the integral vanishes. Thus the contribution along the cut directly related 
to the discontinuity is minus the contribution from the contour integral along the contour 
C which stays far away from the resonance region. 

The distance between the right end of the physical cut and the left end of the other cut 
is 2 while the maximum length of the physical cut is 1 for ?/c = and r = 0. Therefore we 
can reliably calculate the contour integral along C without encountering the other cut. T 
can be expanded in powers of and 1/mi,. The leading term in l/mt is given by 

where Q = rribV — k is the momentum of the s quark and Sf{x) denotes the s quark 



propagator. In momentum space the leading term of Eq. (|2.15| ) in as is proportional to 



(2.16) 



(Q + A;,)2-m2' 
where kb is the residual momentum of the by field of order A. 

For <^ <^ ml it is sufficient to keep the terms at leading order in l/rub only and 
expand Eq. ( p.l5| ) in powers of kb/\Q\. This generates the coefficient functions of the local 



operators in which kb is replaced by the covariant derivatives acting on the fields. The 
leading term independent of kb is given by 



T = 16m^(l + r)^-^, (2.17) 



which reproduces the parton- model result. Here we have used the normalization 
{BW'~f^by\B) = 2mbV^. Eq. ( p.l7|) can be expressed in terms of y as 
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To = -4m^(l + r) (2.18) 

y — 1 + r 

The imaginary part gives rise to a delta function at the endpoint 

W = 21mT = 87rm^(l + r)y^6{l -y-r). (2.19) 



We finally obtain 



^=ro<5(l-2/-r), (2.20) 
dy 

where 

To = W.ncr(m,)r(l + r)(l - r)^ (2.21) 
is the total decay rate in the parton model. 

III. INSTANTON CONTRIBUTION 

We now compute the instanton effect to the decay rate using the dilute gas approxima- 
tion. We emphasize that we have not attempted to calculate other nonperturbative effects 
such as multi-instantons or renormalons. It is possible that they give bigger contributions 
than what we compute. Neverthless our estimate may still characterize a typical size of the 
nonperturbative effects. 

In estimating the contribution we analytically continue the amplitude into the Euclidean 
region of the kinematic variables where is large enough to use the OPE reliably. In the 
background of an instanton (+ anti-instanton) of size p and instanton orientation U located 
at the origin, the Euclidean fermion propagator may be expanded in small fermion mass 



as 
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S±{x, y; p±; U±) = -— ?/^o(a;)^o(y) + S^^\x,y; p±; U±) 

m 

+ mj d'^wS^l\x,w; p±;U±)S^l\w,y; p±;U±) 

+ 0(m2), (3.1) 



where ± denotes instanton, anti-instanton. ipo is the fermion zero mode eigenfunction and 
S±:^ = J2e>o E±{x)'^^E±iy) is the Green function of fermion nonzero modes. 

In evaluating the forward Compton scattering amphtude T, Eq. (|2.15| ), we use the prop- 



agator in Eq. instead of the free propagator. T can be written as 

T = To + Ti„,t, (3.2) 

where the first term is the parton-model amphtude. The second term is the amphtude due 
to instantons of aU orientations U, position z and size p. After averaging over instanton 
orientations Tjnst is given by 



X (fi|6,(x)[A7l(l + 75){5.(X,r;pJ-5o(A)}[/;,7^](l-75)6,(y)|i?), (3.3) 



where D{p) is the instanton density, A = x — y,X = x — Za and Y = y — Za- In Eq. ( p.3|) 
S±{X, Y; p±) is the fermion propagator averaged over instanton (anti-instanton) orientations 
centered at z and 5*0 is the free fermion propagator. Using the MS scheme with ri/ flavors 
of light fermions, D{p) is given by 

45— 5n j 

D{p)=KAHpArf 0";^)"^' (3.4) 

where 



X ^ exp[i - a(l) + 2(n; - l)a(i)] (3.5) 



in which the f3 function at two loops and the running mass at one loop are used and rhi 
are the renormalization-invariant quark masses. In Eq. ( |3.5| ) a{l) = 0.443307 and a (1/2) = 
0.145873. From now on we replace the logarithmic term in D{p) by its value for p = 1/\Q\- 
Corrections to this replacement are negligible since they are logarithmically suppressed. 

In Eq. ( ^.31 ), due to the chirality structure, the nonzero- mode Green function is the 
leading contribution in the 0{ml) and Olms)"^ terms while the zero-mode Green function 
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contributes to 0{mbms) term. However the 0{mbms) term is proportional to k"^. This 
vanishes since the outgoing photon is on the mass shell. Therefore the leading instanton 
contribution comes only from the nonzero-mode Green function. For simplicity we neglect 
the ml /ml subleading term. 

The difference of the single-instanton propagator and the free propagator in momentum 
space after integrating over the instanton orientation, position and size is given by lTl|] 

where 

olO 

iVinst = K—7^'T\Q)k^\ (3.7) 

The exponent of the logarithm is almost unity and we replace it by one. We expect that 



this does not change the result much. Following the same calculation as in Ref. |TT[, we find 
the instanton contribution to T is 

k ■ Ok ■ V O"^ 

T,nst{v ■ k) = 16mj|iVi„3t ^Q^y In ^. (3.8) 

For simplicity we put r = in Eq. ( |3.^ ). The inclusion of nonzero r is straightforward. In 
terms of the rescaled kinematic variables we find 

4iVinst y2 2 

^ins. = ^^^ln[^(l-y)]. (3.9) 

Instanton effects on the B decay depend on the momentum cutoff of the final- 

state photon. The cutoff to remove soft photons should be introduced in experiments since 
it is difficult to isolate a soft photon from the background coming from the subsequent decay 
of, say, K* ^ K + 'j. Experimentally the photon spectrum from B Xg'y is concentrated in 
the region 2.2 GeV < < 2.5 GeV. In CLEO the lower cut is taken at ~ 2.2GeV §. 
Thus we restrict the region of contour integral to yc < y < 1 where yc is the experimental 
cutoff for the scaled photon energy. Then the instanton decay rate is 



Tinstd/c) = To 



inst 



60mP 



10 36 45 20 



y,r (i-yc)' (1-yc)^ (1-1/c 
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(3.10) 



For Uf = 3, the numerical value of K is approximately 126 murhdrhs/A^. Using the 
renormalizat ion-group invariant quark masses m„ = 8.2 ± 1.5 Mev, rhd = 14:A ± 1.5 MeV 
and rhs = 288 ± 48 MeV [0, Eq. {^J§ is written as 



TinstiVc) _ ^ 6.67GeV y f Ay ^ 10 + ^0 



ro \ nib J yrribJ 



. (3.11) 



;i-y,)^ {i-y,r (1-2/er (i-i/c 

The ratio rinst(yc)/ro is shown in Fig. 2 as a function of the cutoff i/c for different values 
of A = 350,400,450 MeV with mb = 5 GeV. For ^ 0.82 that CLEO has chosen, the 
instanton correction is (1.37 ~ 13.1) x 10~^ as A varies from 350 MeV to 450 MeV. Therefore 
the instanton correction is negligibly small compared to other corrections. Note that the 
instanton contribution is much smaller as the cutoff yc decreases. On the other hand, as 
Fig. 2 shows, the contribution is appreciable at yc ^ 0.92. In this region we need a smearing 
to make the perturbation theory valid. However this cutoff is too large to be significant 
experimentally since there is a very small fraction of the rate in this energy window. 

The a^-correction to the decay rate also becomes large as yc approaches 1. Ali and Grueb 
]TI| have examined the correction in detail and have concluded that the a^-correction has to 
be exponentiated near the endpoint of the photon spectrum. Numerically they have found 
that the exponentiation is necessary for yc ^ 0.85. Combined with their result our analysis 
indicates that there is a large theoretical uncertainty for the cutoff yc ^ 0.88. Therefore in 
order to compare experiments with theory in a model-independent way the cutoff yc has to 
be chosen below 0.88. 

We also obtain the instanton contribution to the differential decay rate as 

'^rinst -p A^inst U /'o io\ 

^0 12 n ^-^-^^^ 



Numerically we find 



dy mf (1 -y)' 



1 ^ri„3t ^ / 26.1GeV \3 A ^3 

To dy \ nib ) ^mj {I - yf' ^ " ' 



The ratio in Eq. (|3.13D is shown in Fig.3 for different values of A = 350,400,450 GeV with 



nib = 5 GeV. The ratio also increases sharply near the endpoint of the photon spectrum. 
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Above y ^ 0.87 the instanton correction is appreciable and necessitates a smearing in this 
region. 



IV. DISCUSSION ON KINEMATICS 

We have found that the instanton effect is small in the total decay rate of i? ^ 
as long as Uc is small enough. This is in contrast to the case of i? — > XuCVf. decay. This 
apparently big difference can be understood from kinematics. 

As explained in Section 2 we have obtained the averaged total decay rate by deforming 
the integration contour (from C to C in Fig.l). If the deformed contour is chosen suffi- 
ciently far away from the resonance region we can calculate the decay rate reliably using 
the perturbation theory. It depends on the details of the kinematics whether such contour 
deformation is possible. 

In 5 decay the radius of the deformed contour C is fixed at 1 — ?/c- This radius 

represents the off-shell invariant mass-squared of the final hadrons. Note that this radius 
is independent of kinematic variables. As long as 1 — ?/c ^ Is? /ml the averaged decay rate 
can be calculated reliably at the scale ml{l — y^. The instanton contribution is evaluated 
at this scale. As calculated in Section 3, the contribution is (9(10~^). 

In _B — > XueVf. decay the radius of the deformed contour is given by ^ = {y — q'^){l — y)/y 
where y = 2Ef./'mi, and = [kg + kyY /ml [0. The off-shell invariant mass-squared of 



the final hadrons is equal to m^z. As long as z ^ h?' /m^ we can calculate the averaged 
decay rate using perturbation theory. However unlike the B Xg'-y case the radius z 
depends on how the momentum transferred to the leptonic system is distributed between 
the electron and the anti-neutrino. In particular z vanishes at the boundaries of the phase 
space. Therefore the instanton contribution evaluated at the scale m^z grows rapidly near 
the boundaries y = 1 and y = (f- This means that it is impossible to avoid the resonance 
region unless we introduce a model-dependent cut near the boundaries of the phase space. 
Note that the dependence of the off-shell invariant mass-squared of the final hadrons on 
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the kinematic variables arises only when there are more than one "nonhadronic" particle 
in the final state. In this respect it is interesting to compare our result with the inclusive 



hadronic r decay t ^ + X [0] . In this case the neutrino plays the same role as the 
photon in the B Xs7 decay as far as the kinematics is concerned. The maximum invariant 
mass-squared s = of the final hadrons is fixed at independent of the neutrino energy. 
Since this is much larger than the total inclusive decay rate can be calculated reliably 
using the deformed contour. While the specific kinematic variables under consideration are 
different, the analytic structure and the fact that the maximal invariant mass-squared of the 
final hadrons is independent of other kinematic variables are similar in both B — > and 
inclusive hadronic r decays. 



Indeed the instanton contributions are small in both cases. Nason and Porrati have 
obtained the instanton contribution to the ratio of the hadronic to the leptonic width R^- 
which is given by 

i?;-* / 3. 64A^ 9 m^m^m, ^^^^ 



\ rrir J ml 

Our result for the instanton contribution to the decay rate oi B ^ X^'y is written as 

To \ mb J ml 

for Uc = 0. The fact that both have the same mass dependence and similar numerical factors 
confirms our expectation. 

By the same argument we expect the instanton contribution to the B X^e^e" decay 
rate is similar to that of the B XuCVf. because kinematics and the analytic structure of 
the forward Compton scattering amplitudes are the same. 

V. CONCLUSION 

In this paper we have calculated the instanton effect on the decay rate oi B ^ ^s7. 
Unlike B X„ez7e decay we have found that the instanton correction is negligibly small. 
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We have also observed that the kinematics and the analytic structure in this case is quite 
similar to the inclusive hadronic r decay, in which the instanton correction is also small. 

It is interesting to note that the contribution of nonperturbative effects like instantons 
depends on kinematics. It depends on whether kinematics allows to choose a contour in the 
complex plane of the relevant variables to avoid the resonance region and to use perturbation 
theory reliably. In the inclusive semileptonic B decay it is impossible to choose a deformed 
contour to avoid the resonance region. Therefore we need a smearing prescription which 
limits the theoretical prediction in a model-independent way. 

When we probe new physics from the radiative inclusive B decay, it is therefore suffi- 
cient to examine only those corrections coming from the hadronic matrix elements and the 
radiative QCD corrections. The latter part is ~ 30% of the parton-model result. On the 
other hand since the QCD coupling constant is sensitive to the choice of the renormalization 
scale fi at one-loop order, it is necessary to set the scale /i precisely. This requires QCD 
corrections beyond one-loop. Partial attempt along this direction has been made recently 
by Bur as et.al. |TH . 
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FIGURES 

FIG. 1. Anaytic structure of T in the complex y plane. The point P at y = t/c is the lower 

cutoff of the photon energy. The contour integral along C is minus the contour integral along C. 

FIG. 2. The ratio rinst(yc)/ro for A = 350,400,450 MeV with = 5 GeV. 
FIG. 3. The ratio (l/ro)drinst/d2/ for A = 350,400,450 MeV with = 5 GeV. 
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